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0.  ABSTRACT 

Properties  of  the  Shevhart  X-chart  for  controlling  the  mean  of  a 

process  with  a  normal  distribution  are  investigated  for  the  situation  where 
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the  process  variance  a  must  be  estimated  from  initial  sample  data.  The 

—  2 

control  limits  of  the  X-chart  depend  on  the  estimate  of  a  and  thus ,  unlike 

the  case  when  is  known,  the  X-chart  is  not  equivalent  to  a  sequence 
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of  independent  tests.  When  a  is  estimated  the  distribution  of  the  run 

length  is  not  geometric  and  cannot  be  characterized  simply  in  terms  of  the 

probability  of  a  signal  at  a  given  point.  The  average  run  length  (ARL)  for 

the  X-chart  is  expressed  in  terms  of  an  integral  involving  the  normal  cdf, 

and  it  is  shown  that  the  chart  signals  with  probability  one,  but  the  ARL 

2 

may  not  be  finite  if  the  size  of  the  sample  used  to  estimate  a  is  suffi¬ 
ciently  small.  In  addition,  certain  bounds  for  the  ARL  are  also  derived. 

Numerical  integration  is  used  to  show  that  the  effect  of  using  small  sam- 

2 

pie  sizes  in  estimating  o  is  to  increase  the  ARL  and  the  variance  of  the 
run  length  distribution. 

*This  work  was  supported  in  part  by  Grant  No.  DAAG-29-78-G-0172  from  the 
U.S.  Army  Research  Office. 

KEY  WORDS:  X-Chart;  Shewhart  Chart;  Process  Control;  Control  Chart; 
Estimated  Variance;  Economic  Model;  Average  Run  Length. 
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1.  INTRODUCTION 

Control  charts  have  long  been  used  to  detect  changes  in  the  distri¬ 
bution  of  observations  taken  from  a  process  that  operates  continuously 
over  time.  Various  criteria  have  been  proposed  to  measure  the  perfor¬ 
mance  of  a  control  chart  when  it  is  used  to  control  a  process.  These 
criteria  aid  in  the  design  of  charts  to  meet  specific  objectives  and  also 
serve  as  a  basis  for  comparing  several  competing  procedures.  The  purpose 
of  this  paper  is  to  examine  some  of  these  criteria  in  terms  of  their  ap¬ 
propriateness  with  special  reference  to  the  Shewhart  X-chart  for  control¬ 
ling  the  process  mean  when  the  process  variance  in  unknown. 

Consider  the  situation  where  an  independent  sequence  X^,  Xj,  •  •• 
of  random  samples  is  taken  at  regular  intervals  from  the  process, 

■  (Xii*  xi2’  * ' *  ,Xin^  ^eing  the  sample  of  n  observations  taken  on 
the  1th  occasion.  In  most  cases  Che  distribution  of  the  observations  will 
depend  on  one  or  more  parameters  represented  by  8,  which  may  be  vector- 
valued.  The  basic  problem  is  to  detect  quickly  any  change  in  the  value 
of  9  by  using  the  values  of  all  observations  taken  up  to  the  current  time. 
In  some  situations  the  Initial  control  value,  say  9q,  for  9  will  be  known 
or  specified  in  advance.  In  other  situations  one  or  more  components  of 
9g  may  be  unknown,  and  it  is  necessary  to  estimate  these  values  using  the 
observed  data. 

A  process  control  procedure  for  detecting  shifts  away  from  9q  is  a 
procedure  that  at  time  i  decides  on  the  basis  of  X^,  Xj ,  ...  ,X^  whether 
to  continue  sampling  or  to  stop  and  signal  that  a  change  has  occurred. 

As  long  as  the  value  of  9  remains  constant  at  9g  we  want  to  continue 
sampling,  but  as  soon  as  a  significant  shift  occurs  we  want  to  stop. 

If  N  represents  the  random  time  at  which  the  procedure  stops  and  signals, 
then  it  is  the  distribution  of  N  that  determines  the  properties  of  the 
procedure.  In  process  control  terminology,  N  is  called  the  run  length.. 


As  an  example  of  a  process  control  procedure,  consider  the  standard 
Shewhart  X-chart  for  detecting  shifts  in  the  mean  u  of  a  process  having 
a  normal  distribution.  This  procedure  signals  that  a  shift  has  occurred 
at  the  first  value  of  i  for  which 

h  s  *0  + 

where  X^  is  the  mean  of  X^,  X^*  •••  »Xin»  wq  is  the  target  or  control 
value  for  the  process  mean,  k  is  a  constant,  and  a  is  the  process  standard 
deviation.  In  what  follows  it  is  convenient  to  express  the  process  mean 
in  terms  of  $  ■  (y  -  y^'/n/o,  the  case  y  -  Uq  then  corresponding  to  6  *  0.  • 
If  both  yQ  and  o  are  known  (or  specified) ,  then  the  distribution  of  the 
run  length  is  geometric  with  parameter 

p(<5)  -  P(Xi  s  uQ  +  ko)  -  1  -  *(«  +  k/n) 

since  we  are,  in  effect,  performing  a  sequence  of  "independent  tests"  where 
the  ith  test  involves  only  the  ith  sample.  Here  $(•)  is  the  cumulative  dis 
trlbution  function  of  the  standard  normal  distribution.  The  parameter  p(6) 
completely  determines  the  run-length  distribution  and  thus  the  behavior  of 
this  procedure  can  be  summarized  conveniently  and  naturally  in  terms  of 
this  probability  of  a  signal  on  a  given  occasion.  In  addition,  if  we  let 
^  b*  a  value  of  <5  that  represents  a  shift  of  interest,  then  the  quantities 
p(0)  and  l-p(6^)  can  be  thought  of  as  the  probabilities  of  errors  of  the 
first  and  second  type,  respectively,  for  testing  the  statistical  hypothesis 
H:  5  •  0  against  A:  $  ■  <5^.  This  formulation  using  the  concepts  and  ter¬ 
minology  of  hypothesis  testing  is  very  natural  to  statisticians. 

The  drawback  underlying  p(<5)  as  a  general  parameter  for  process 
control  procedures  for  detecting  changes  in  5  is  that  it  may  not  be  mean¬ 
ingful  for  more  complicated  procedures  such  as  cumulative  sum  (CUSUM)  con¬ 
trol  charts,  moving  average  charts,  or  Shewhart  X-charts  with  estimated 
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variance.  For  example,  the  probability  that  a  CUSUM  chart  signals  on  oc¬ 
casion  i  is  not  independent  of  i  since,  for  this  procedure,  the  decision 
at  stage  i  depends  on  X^,  X^,  ...  ,X^_^  in  addition  to  X^.  When  the  pro¬ 
cess  variance  is  unknown  and  must  be  estimated  from  preliminary  samples, 
the  procedure  no  longer  entails  a  sequence  of  Independent  tests  and  the 
run-length  distribution  cannot  be  characterized  by  the  probability  of  a 
signal  at  a  given  point.  Thus,  for  general  process  control  procedures, 
some  additional  measures  or  parameters  are  needed  to  describe  the  run- 
length  distribution.  One  commonly  used  parameter  is  the  mean  E.(N),  called 
the  average  run  length  (ARL).  Another  useful  parameter  is  P  (N  s  t) ,  where 
t  is  a  fixed  value. 

One  objective  of  this  paper  is  to  determine  which  of  the  measures 
associated  with  the  run-length  distribution  are  most  important  in  deter¬ 
mining  the  characteristics  of  the  procedure.  This  is  done  through  the 
development  of  a  simple  economic  model  in  Section  2.  Another  objective 
is  to  examine  the  properties  of  the  X-chart  when  the  variance  of  the  pro¬ 
cess  must  be  estimated.  In  Section  3,  various  methods  of  estimating  the 
variance  are  discussed  and  the  effect  of  estimated  variance  on  the  ARL  is 
investigated.  Properties  of  the  X-chart  are  usually  calculated  under  the 
assumption  that  a  is  known,  but  our  results  indicate  that  the  actual 
properties  of  the  X-chart  when  o  is  estimated  can  be  quite  different  from 
the  properties  for  the  case  when  a  is  known.  This  distinction  is  often 
Ignored  in  the  literature.  Section  4  extends  the  results  of  Section  3  to 
the  two-sided  X-chart. 

2.  AN  ECONOMIC  MODEL 

To  determine  the  characteristics  of  the  run  length  distribution  that 
are  important  determinants  of  the  properties  of  a  control  chart,  we  may  use 


an  economic  model  that  evaluates  the  Impact  on  Income  and  costs  associated 
with  the  use  of  a  particular  control  chart  (see  Duncan  (1956)).  The  economic 
models  that  have  been  developed  In  the  literature  usually  express  the  long¬ 
term  average  net  income  per  unit  time  for  the  process  as  a  function  of  vari¬ 
ous  process  parameters  and  of  various  parameters  associated  with  the  control 
procedure.  The  usual  objective  in  the  design  of  control  procedures  is  to 
maximize  this  net  income. 

A  simplified  economic  model  will  now  be  developed  to  show  how  the  character¬ 
istics  of  the  run-length  distribution  enter  into  the  model.  This  model  is 
for  a  general  process  control  procedure  for  controlling  the  general  parameter 
6  in  contrast  to  models  in  the  literature  which  are  specific  to  the  procedure 
and  parameter.  Assume  that  the  process  is  operating  in  control  when  the  con¬ 
trol  procedure  is  first  applied  and,  for  simplicity,  let  the  unit  of  time 
correspond  to  the  time  interval  between  samples.  At  some  random  time  Tq 
(counted  from  the  start  of  the  procedure) ,  the  process  mean  shifts  away  from 
the  control  value  0q.  Assume  for  simplicity  that  there  is  only  one  possible 
out-of-control  state  represented  by  the  value  9^.  After  the  shift  from  0q 
to  0^  the  process  remains  at  0^  until  the  control  procedure  detects  the  out- 
of-control  situation  and  rectifying  action  is  taken  to  return  0  to  0q.  Ano¬ 
ther  in-control-out-of-control  production  cycle  then  starts  at  the  instant 
when  9  returns  to  9q.  Thus,  the  long-term  application  of  a  control  procedure 
involves  an  infinite  sequence  of  these  in-control-out-of-control  cycles. 

By  considering  income  lost  due  to  operating  the  process  out  of  control 
as  a  cost,  the  objective  of  maximizing  the  long-term  average  net  income  per 
unit  time  can  be  reformulated  as  a  problem  of  minimizing  the  long-term  aver¬ 
age  cost  per  unit  time.  In  each  production  cycle  the  primary  costs  are  those 
due  to  false  alarms,  to  operating  out  of  control,  and  to  sampling.  The  long¬ 
term  average  cost  per  unit  time  is  obtained  by  finding  the  expected  cost  per 
cycle  and  then  dividing  the  latter  by  the  expected  length  of  a  cycle  (see 
Johns  and  Hiller  (1963)). 
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Let  T  represent  the  time  (counted  from  the  start  of  the  cycle)  when 
the  procedure  detects  Che  shift  and  8  is  returned  to  8q.  Then  T  is  Just 
the  random  cycle  length.  Let  W  represent  the  number  of  false  alarms  dur¬ 
ing  the  part  of  the  cycle  when  the  process  is  in  control.  The  cost  per 
cycle  can  then  be  represented  by 

C  -  Cl(W)  +  c2(T-Tq)  +  c3(T) , 

where  the  functions  c^,  c^,  and  c^  are  given  by 

c^(w)  =»  cost  due  to  w  false  alarms  during  a  cycle, 

c2(t)  ■  cost  due  to  operating  out  of  control  for  t  time  units, 

c^Ct)  *  cost  of  sampling  during  a  cycle  of  length  t. 

The  long  run  average  cost  per  unit  time  is  therefore 

A  -  E(C)/E(T). 


It  seems  reasonable  to  assume  that  the  cost  of  false  alarms,  c^(W),  is 

directly  proportional  to  W  so  that  c^(W)  *  a^W  where  a^  is  a  constant.  In 

this  case  E(c.(W))  is  approximately  equal  to  a1E(Tn)/E  (N) .  This  approxima- 
l  l  U  9q 

tion  can  be  justified  as  follows.  Let  N^,  N2,  ...  be  the  sequence  of  run 
lengths  between  false  signals  during  the  in-control  part  of  the  cycle,  and 
let  ^  be  the  run  length  after  that  would  be  required  for  a  signal  if  the 
shift  from  0q  to  9^  does  not  take  place.  Then  ^  has  the  same  distribution  as 
N^.  Thus 


W 

r  n 

i«l 


i 


<  T„ 


W-l 

SIN 

i«l 


and,  since  E. 

90 


E(W)E.  (N), 


E(Tq)/E9  (N)-l  S  E(W)  <  E(T0)/E@  (N)  . 


In  many  cases  it  may  be  reasonable  to  assume  also  that  C2(t)  Is  of  the 
form  a2t  where  a2  is  a  constant.  Neglecting  the  fact  that  a  shift  can 
occur  between  samples  and  assuming  that  the  expected  time  to  detect  the 
shift  does  not  depend  significantly  on  how  long  the  control  procedure 
has  been  running  in  control,  E^CT-Tq))  is  approximately  equal  to 
a2E0  (N).  The  expected  cycle  length  E(T)  is  approximately  E(TQ)  + 

E0  (N) .  The  sampling  cost  c^OT)  depends  only  on  T  and  n.  The  long¬ 
term  average  cost  per  unit  time  then  can  be  expressed  as 


a  Ef T  WE  +  a  E  fKrt  ■+•  Err-  fTl 


E(T0)  +  Eg  (N) 


Thus  A  depends  on  the  run-length  distribution  primarily  through  the  ARL 

values  E  (N)  and  E.  (N) .  This  suggests  that,  if  the  economic  model  is 
90  91 

a  reasonable  representation  of  the  actual  cost  and  loss  structure,  the 
primary  measure  of  the  effectiveness  of  control  procedures  is  the  ARL. 

There  may  be  situations,  however,  where  the  function  C2(t)  is  not 
linear  in  t.  For  example,  a  small  amount  of  low  quality  material  pro¬ 
duced  while  operating  out  of  control  may  not  be  particularly  damaging, 
but  if  the  out-of-control  time  exceeds  a  certain  value  then  a  complete 
production  run  may  be  jeopardized.  In  this  case  the  cost  function  could 

take  the  form  »  0  if  t  s  t^  and  CjCt)  *  a^  if  t  >  t^.  In  other 

2 

cases  a  cost  function  of  the  form  C2(t)  »  t  might  be  appropriate.  Thus, 

2 

in  addition  to  the  ARL,  parameters  such  as  P  (N>t.)  or  E  (N  )  may  be 

61  1  91 

of  interest  in  particular  applications. 

3.  ONE-SIDED  X- CHARTS  WITH  UNKNOWN  VARIANCE 

In  this  section  we  consider  the  one-sided  Shewhart  X-chart  for  de¬ 
tecting  shifts  in  the  positive  direction  in  the  mean  of  a  process  having 
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a  normal  distribution.  A  one-sided  X-chart  for  detecting  shifts  in  the 
negative  direction  can  be  obtained  by  using  a  lower  limit  and  the  results 
obtained  in  this  section  can  be  applied  by  replacing  6  with  -6  where,  as 
before,  5  ■  (ji  -  Uq)  */n 7a.  The  two-sided  X-chart  will  be  discussed  in  Sec¬ 
tion  4. 

When  the  variance  of  the  process  is  unknown,  some  fora  of  estimate  of 
the  variance  must  be  used  in  the  construction  of  the  control  limits  for  the 
chart.  The  estimate  of  the  variance  can  be  obtained  from  the  information 
in  each  sample,  or  from  an  initial  sequence  of  samples,  or  from  a  combina¬ 
tion  of  the  two  approaches  where  the  initial  variance  estimate  is  updated 
as  each  new  sample  is  taken.  The  variance  estimates  are  traditionally  based 
on  either  the  sample  variance  or  the  sample  range.  Since  the  sample  vari¬ 
ance  is  more  efficient  than  the  sample  range,  we  will  restrict  attention 

to  situations  where  sample  variances  are  used.  We  now  look  at  the  impli- 

2 

cations  of  the  various  ways  of  estimating  a  as  they  bear  on  the  properties 
of  the  X-chart. 

If  the  sample  variance 

2  n  —  2 
s  -  Z  (X  -X  )  / (n-1) 

for  stage  i  is  used  as  the  variance  estimator  at  stage  i,  then  the  procedure 
that  signals  at  the  first  i  for  which 

Xf  >  u0  +  ksi 

is  a  sequence  of  independent  tests.  The  parameter 

p*(5)  -  P(X±  >  uQ  +  kSi) 
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can  be  calculated  using  the  noncentral  t-distribution  for  any  value  of  6 

and  the  ARL  of  the  procedure  is  1/p* (6).  This  procedure  is  not  used  widely 

in  practice  because  n  is  usually  too  small  to  provide  a  reasonable  estimate 
2 

of  a  when  each  sample  is  used  independently.  In  addition,  the  procedure 
must  be  carried  out  in  its  standardized  version 

<Xt  -  p0)/S1  >  k 

in  order  to  have  the  same  upper  limit  at  each  stage,  and  this  may  inhibit 
a  visual  interpretation  of  the  plots  since  the  chart  is  not  showing  the 
actual  sample  means. 

2 

A  second  approach  to  the  estimation  of  o  involves  the  use  of  the  same 
initial  estimate  for  each  stage.  In  most  cases  this  initial  estimate  is 
based  on  a  sequence  of  r  initial  samples,  each  of  size  m.  In  addition  to 
providing  a  variance  estimate,  these  r  samples  can  be  examined  further  to 
determine  whether  the  process  was  in  control  with  respect  to  variance  at 
the  time  that  these  samples  were  taken.  If  the  r  sample  variances  are 
pooled  to  obtain  the  estimator 

2  r  2 
S  -  Z  S“7r 

1-1 

then  the  procedure  signals  at  the  first  i  for  which 

Xt  >  m0  +  kS. 

This  approach  has  the  advantage  that  the  chart  shows  the  actual  sample  means. 
Since  the  same  S  is  used  at  each  point,  the  procedure  no  longer  entails  a 
sequence  of  independent  tests,  and  the  run-length  distribution  will  not  be 
geometric.  In  this  case  the  parameter 

p+  (5)  -  PO^  >  uQ  +  kS) 

will  not  suffice  to  characterize  the  run-length  distribution. 
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2 

As  an  alternative  to  the  pooled  estimate  of  a  from  a  sequence  of  ini- 

2 

tial  samples,  one  large  initial  sample  can  be  taken  to  estimate  a  .  A  dis¬ 
advantage  of  this  approach  is  that  it  is  not  possible  to  test  for  control 

of  the  variance  at  the  time  the  sample  was  taken. 

2 

A  third  approach  to  estimate  o  uses  some  kind  of  initial  estimate  and 
then  updates  the  estimate  as  each  new  sample  is  taken.  Although  this  method 
appears  reasonable  as  long  as  the  process  variance  is  fixed,  additional  com¬ 
putational  effort  would  be  required  and  the  control  limits  would  change  after 
each  sample.  In  addition,  the  properties  of  such  a  procedure  are  quite  dif¬ 
ficult  to  determine. 

A  common  practice  in  setting  up  a  control  chart  is  to  take  k  ■  3/t^T  re- 
2 

gardless  of  whether  a  is  known  or  estimated  from  a  small  sample  (for  example, 

see  Duncan  (1974)).  Various  authors  (cf.  King  (1954),  Proschan  and  Savage 

(1960),  Hillier  (1964),  and  Yang  and  Hillier  (1970))  have  considered  the 

2 

problem  of  choosing  k  when  a  is  estimated.  Yang  and  Hillier  (1970)  proposed 

+"  2 

that  k  be  chosen  to  give  a  specified  value  of  p  (0) .  If  the  estimator  S  has 
v  *  r(m-l)  degrees  of  freedom,  then  the  choice 

k  -  t.  (v)//n 
1-a 

gives  p+  (0)  =  a,  where  t-  (v)  is  the  1-a  fractile  of  the  t-distribution 

l-a 

having  v  degrees  of  freedom.  By  choosing  t,  (v)  corresponding  to  a  *  .00135, 

1-a 

the  value  of  p  (0)  will  be  the  same  as  the  value  of  p(0)  for  k  *  3/^n.  Of 
course,  t^_^(v)  can  be  chosen  to  give  any  other  desired  value  of  p*(0) .  The 
problem  with  this  approach  is  that  the  run-length  distribution  and  ARL  are  still 
unknown.  It  will  be  shown  later  in  this  paper  that  the  actual  ARL  for  the  pro¬ 
cedure  using  S  as  an  estimator  for  o  can  differ  substantially  from  the  ARL 
calculated  under  the  assumption  that  p  (0)  »  p(5). 
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2  2 

If  the  estimator  S  for  a  is  a  sample  variance  or  a  pooled  sample  var- 

2  2 

iance  based  on  v  degrees  of  freedom,  then  \>S  /a  has  a  chi-squared  distribu¬ 
tion  with  v  degrees  of  freedom  when  the  process  distribution  is  normal.  Con¬ 
sequently  the  run  length  N  satisfies 

P (N  >  t)  -  P(X  <  uQ  +  kS,  i  -  1,  2,  ...  ,t) 

00 

«  /  (-5  +  —  /a)  f  (s)ds 

0  °  S 

*  ry-r  /  $t(-5  +  c/T)uS^V  2^e  Udu  (3.1) 

r(hv)  0 

where  f  (s)  is  the  density  of  S,  and  c  *  k/2n/v.  Unfortunately  the  integral 

O 

in  (3.1)  cannot  be  expressed  in  a  simpler  closed  form.  Nevertheless,  certain 
results  for  the  run-length  distribution  can  be  established  using  (3.1).  These 
results  are  given  in  the  following  theorems,  and  the  proofs  are  given  in  the 
Appendix. 

THEOREM  1.  For  all  values  of  k  and  5,  P(N  <  “)  *  1  and  lim  P(N  >  t)  =>  (-<$+k/^T) 

v-w 

This  theorem  implies  that  the  X-chart  using  an  initial  estimate  of  a  eventually 

will  signal,  whatever  the  values  of  u  and  a  may  be.  Moreover,  as  the  number 

2 

of  degrees  of  freedom  for  the  initial  estimate  of  a  grows  larger,  the  distri¬ 
bution  of  the  run  length  converges  to  the  geometric  distribution  for  the  case 
2 

where  a  is  known.  The  results  given  in  this  theorem  hold  without  the  assump¬ 
tion  of  normality. 

The  ARL  of  the  X-chart  can  be  expressed  as 

00 

E(N)  -  Z  P(N  >  t). 
t-0 

Using  this  expression  the  following  theorem  can  be  proved. 
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THEOREM  2. 

(a)  If  c  <  Jl  then  E(N)  <  «  for  all  5. 

(b)  If  c  >  J2  then  E(N)  »  «  for  all  5. 

(c)  If  c  ■  J2  then  E(N)  <  »  for  6  >  0  and  E(N)  *  «  for  S  s  0. 

(d)  If  E(N)  <  «  then 


E(N) 


1 

r(hv) 


®  4(v-2)  -u 

/  — - 2 - -=r 

0  1  ~  $(-5  +  c/u) 


du. 


(e)  lim  E(N)  »  Cl  -  $(-6  +  k^)]“L. 


(3.2) 


This  theorem  gives  the  rather  surprising  result  that  under  certain  conditions 

the  ARL  of  the  X-chart  with  estimated  variance  is  not  finite.  The  condition 

2 

c  >  Jl  is  equivalent  to  the  condition  v  <  k  n.  Thus  if  v,  the  number  of  de- 

2 

grees  of  freedom  in  the  estimate  of  a  ,  is  sufficiently  small,  the  ARL  is  not 

finite.  In  particular,  if  k  is  taken  as  the  3/ Sn  (the  standard  3-sigma  limits) 

then  c  >  /?  if  v  <  9.  If  we  take  k  *  t,  (v)//ii,  following  the  suggestion  in 

1-ct 

2 

Yang  and  Hillier  (1970), then  c  >  /2  whenever  v  <  t^_a(v)  (e.g.  when  a  *  0.001 
the  ARL  is  not  finite  if  v  S  14).  However,  if  v  is  large,  say  100,  then 

c  <  1/2"  if  k/n  <  9.95  and  the  ARL  will  be  finite  for  values  of  k  and  n  that 

are  likely  to  be  used  in  applications. 

The  following  general  bound  for  the  run-length  distribution  can  be  obtained. 
THEOREM  3.  For  all  values  of  5,  k,  and  v, 

P(H  S  t)  S  1  -  (l-p+(5))C  . 

‘Thus  the  run-length  distribution  is  stochastically  larger  than  a  geometric 
distribution  with  parameter  p+($).  This  theorem  holds  regardless  of  the  form  of 


J 
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the  underlying  distribution  of  the  observations,  although  computation  of  p  (5) 
may  be  difficult  for  non-normal  distributions.  If  the  underlying  distribution 
is  normal  then 

pt(6)  -  P(t(v,6)  >  k/n)  , 

where  t(v,S)  has  a  noncentral  t-distribution  with  \>  degrees  of  freedom  and 
noncentrality  parameter  5.  Suppose  k  is  chosen  as  k  *  t^  a(v)/^n,  following 
the  suggestion  of  Yang  and  Hillier  (1960).  Then,  when  the  process  is  in  con¬ 
trol,  the  geometric  distribution  with  parameter  a  serves  as  a  stochastic 
lower  bound  for  the  run-length  distribution.  The  theorem  also  provides  a 
lower  bound  for  the  ARL  since 

E(N)  >  l/p+(p)  . 

Sharper  bounds  for  the  ARL  under  the  assumption  of  normality  are  given 
in  the  following  theorems. 

THEOREM  4.  If  0  <  c  <  Si  and  6*0,  then 


THEOREM  5.  If  0  <  c  <  /I  and  6  >  0,  then 


4 


U(<5)  r1  P(X2(v)  5  2<52/c2)  +  D  max  [A-fiB,  lSA+((2/Tr)ls-is5)B] 

S  E(N)  £  2P(X2(v)  S  252/c2)  +  D{A+(1-S)B},  (3.4) 

where 

D-  (^)Jsr(Ji(v+l))(l-lSc2)'is(v+1)exp{52/(2-c2)}  . 

V  -  2  u  5 

A  -  vc  Z  a,P(x  (v+l-j)>d),  B  -  2(1-Jjc  )*  Z  b  P(x  (v~j)>d), 
j-0  J  j-0  J 

2  2  2  -1 
d  -  2«*{c  (l-%c  )  )  1  , 

aj  *  (-l)j  {Scd-ljc2)"54^  (j!r(iS(v+2-j))}~1  , 
b  -  (-l)j  {Sc(l-JscW  {j  !r(%(v+l-j))  r1 

2 

and  where  x  (v)  represents  a  chi-square  random  variable  having  v  degrees  of 
freedom. 

THEOREM  6.  If  c  -  JE  and  6  >  0,  then 

A^/$($)  +  K  max[vA3-62A2,  ^vA^  +  (Ojir)**  -  ^SjA^] 

£  E(N)  £  2A1  +  K[vA3  +  6(1-6)A2]  ,  (3.5) 

where 

K  -  ^  -  P(x2(v)«2), 

A2  -  P(x2(2v)  <  262),  A3  -  P(x2(2v  +  2)  >  2«$2)  . 

Thus  far  we  have  examined  the  run-length  distribution  under  the  assump¬ 
tion  that  6  is  constant  for  each  sample.  In  many  situations  the  shift  away 
from  the  control  value  may  take  place  gradually  as  a  drift  over  a  period  of 
time.  Suppose  that  5^  is  the  value  of  5  when  the  sample  is  taken  so  that 
5  *  (5^,  i$2,  ...)  represents  the  set  of  values  of  6  corresponding  to  the 


sampling  occasions.  We  now  Investigate  the  distribution  of  N  as  a  function 
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THEOREM  7.  If  <5  and  6  *  are  two  sequences  such  that  6.  <  6!  for  all  1,  then 

—  -v  11 

P5(N  <;  t)  S  PS,(N  <;  t). 

The  point  of  the  theorem  is  that,  whatever  the  form  of  the  deviation  from 
the  control  value,  the  larger  the  shift  the  smaller  the  detection  time.  The 
computation  of  the  distribution  of  N  for  arbitrary  6  is  difficult.  However, 
if  6  is  bounded  as  5  s  5i  s  6,  for  all  i,  then  letting  5  *  C 6_,  5_,  ...)  and 
o*  (6,  5,  ...)  we  have  that 

PS(N  £  t)  s  P  (H  t)  5  P^(N  s  t). 

Thus  the  run-length  distribution  under  the  assumption  of  a  gradual  drift  can 
be  bounded  by  run-length  distributions  for  constant  process  mean.  The  con¬ 
clusion  of  this  theorem  can  be  shown  to  hold  under  the  general  condition  that 
the  observations  are  from  a  location  parameter  family  with  location  parameter 
6. 

In  order  to  determine  the  effect  of  estimated  variance  on  the  properties 
of  the  X-chart,  P(N  >  t)  and  the  ARL  were  calculated  using  numerical  integra¬ 
tion.  For  the  case  where  n  *  4  and  k  *  3/^n  (the  standard  three-sigma  limits). 
Table  1  gives  the  ARL  and  Table  2  gives  P(N>t)  for  several  values  of  v  and 
6/v'n  *  (u~Uq)/ct.  From  Table  1  we  see  that  the  ARL  is  larger  for  small  values 
of  v  than  it  is  for  the  case  where  a  is  known  (v*®).  From  Table  2  we  see  that 
P(N>t)  is  smaller  than  expected  when  5*0  and  t  *  100  but  larger  than  expected 
for  6  >  0.  For  example,  if  v  ■  29  the  probability  of  stopping  before  100 
observations  when  in  control  is  0.2038  compared  to  0.1264  when  v  *  ®.  But  if 
5/  Jn  ■  1,  the  probability  of  taking  more  than  20  samples  to  detect  the  out-of¬ 
control  situation  is  0.0796  instead  of  0.0316  when  v  *  ®.  One  effect  of  using 
an  estimated  variance  is  to  increase  the  variance  of  N.  This  explains  the 
fact  that,  when  5*0,  E(N)  increases  as  v  decreases  while  F(N  >  100)  decreases 


as  v  decreases. 
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The  ARL  values  are  larger  for  small  values  of  v,  and  Che  effect  of 
using  limits  based  on  the  t-distribution  in  place  of  the  three-sigma  limits 
is  to  increase  the  ARL  values  even  further.  This  suggests  that  if  the 
desired  ARL  values  are  those  calculated  for  three-sigma  limits  when  a  is 
known,  then  the  limits  for  the  case  where  o  is  estimated  should  be  reduced 
rather  Chan  increased.  (Note  that  this  is  contrary  to  the  recommendation 
of  Yang  and  Hillier  (1970)). 

An  additional  complication  arises  in  the  estimation  of  o  when  r  initial 
samples  of  size  m  are  tested  for  control  with  respect  to  variance  before 
being  used  to  estimate  o.  The  usual  practice  is  to  set  up  a  control  limit 
and  discard  any  sample  from  among  the  r  samples  that  exceeds  the  control 
limit.  If  we  let  S*  be  the  resulting  estimate  of  a  then  S*  and  S  have  dif¬ 
ferent  distributions.  The  run  length  N*  for  the  X-chart  in  this  case  satis¬ 
fies 

09 

P(N*  >  t)  -  /  Pt(Xi  S  uQ  +  ks|S*-s)fs*(s)ds. 

Now  S  based  on  the  original  sample  is  stochastically  larger  than  S*  and 
(assuming  k  >  0)  P^  £  y^  +  ks|S*  ■  s)  is  increasing  in  s  so  that  P(N*>  t) 

S  P(N  >  t) .  This  implies  that  E(N*)  £  E(N).  Thus  the  effect  of  testing  for 
control  when  the  Initial  samples  are  taken  is  to  reduce  the  ARL. 
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4.  RESULTS  FOR  TWO-SIDED  SHEWHART  CHARTS 

In  some  situations  it  is  desirable  to  detect  shifts  in  either  direction 
from  the  control  value  Uq.  For  this  problem  it  is  not  necessary  that  the 
control  limits  be  symmetric  about  Uq,  but  symmetric  limits  usually  are  used 
in  practice  and  we  consider  this  case  for  simplicity.  If  an  initial  estimate 
S  based  on  v  degrees  of  freedom  is  used  to  estimate  o,  then  the  two-sided 
X-chart  signals  at  the  first  i  for  which 

!xt  -  ji0I  *  ks. 

Let  N,  N' ,  and  N"  denote  the  run  lengths  of  the  one-sided  chart  for  posi¬ 
tive  deviations,  the  two-sided  chart,  and  the  one-sided  chart  for  negative 
deviations,  respectively.  If  the  three  procedures  are  applied  simultaneously 
using  the  same  value  of  k  then  it  follows  that 

N'  -  mintN.N"}. 

Using  the  results  and  approach  of  Theorem  1  for  the  one-sided  case  we  have 
that  P(N'  <  »)  -  1  and  lim  P(N  <  t)  -  [$(-6  +  k^)  -  *(-«  -  k>£)  ]C. 

Corresponding  to  THEOREM  2  for  the  one-sided  case  we  have  the  following 
result  for  the  two-sided  case. 


THEOREM  8. 

(a)  If  c  <  Jl  then  E(N')  <  «. 

(b)  If  c  >  then  E(N')  ■  •. 

(c)  If  c  -  ^  then  E(N')  <  -  for  5  +  0  and  E(N')  -  »  for  6  -  0. 

(d)  If  E(N')  <  «  then 


The  stochastic  bound  for  the  run-length  distribution  based  on  the 
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geometric  distribution  also  holds  for  the  two-sided  case  when  the  parameter 
P  "*"(<$)  is  replaced  by 

p'(5)  -  P( |X±  -  p0 { /S  2  k). 

When  the  underlying  distribution  is  normal, 

p'(5)  ■  P(jt(v,6)j  ki^u) . 

When  5  »  0,  the  expression  for  E(IT')  given  in  Theorem  8(d)  reduces  to 
E(N)/2  where  E(N)  is  given  by  Theorem  2(d).  Thus  bounds  for  E(N')  for  the 
case  0  <  c  <  /T  and  6*0  can  be  obtained  from  the  bounds  for  E(N)  given  in 
Theorem  4. 

Since  N'  ■  min(N,lT)  it  follows  that  an  upper  bound  for  E(N')  is  given 
by  min(E(N) ,E(N")) .  When  0  <  c  <  ^2  and  5  >  0  an  upper  bound  for  E(N)  is 
given  by  Theorem  5.  When  0  <  c  <  /£"  and  5  <  0  an  upper  bound  for  E(N")  can 
be  obtained  from  Theorem  5  by  using  -5  in  the  bound  for  E(N) . 

Corresponding  to  Theorem  7  for  the  one-sided  case,  we  have  the  following 
theorem  concerning  the  behavior  of  N'  as  a  function  of  5. 

» 

THEOREM  9.  If  5  and  6'  are  two  sequences  such  that  1 6^ [  z  [ 6 ^ [  for  all  i, 
then  P5(N  £  t)  s  P5,(N  s  t) . 

This  result  holds  for  any  underlying  distribution  which  is  symmetric  and 


unimodal . 


APPENDIX 


PROOF  OF  THEOREM  1.  For  any  positive  integer  t 


P(N  >  t)  =  /  $t(-S  +  —  /n)  dF_(s).  (Al) 

„  0  5 


Since  $C(-5  +  ^  i/n)  +  0  as  t  +  *  and  $(-5  +  —  v'n)  £  I  for  every  s,  it  follows 
from  the  dominated  convergence  theorem  that 


lim  P(N  >  t)  -  /  lim  $C(-5  +  —  «/5)dF  (s)  *  0. 
t-**>  0  t-**>  a 

ICS  r— 

Hence  P(N  <  ®)  *  1.  Since  S  o  a.s.  as  \i  +  »  and  $(-6  +  —  /n)  is  contin¬ 


uous  and  bounded,  it  follows  that 


lim  P(N  >  t)  -  lim  I  $C(6  +  —  Vn)dF  (s)  -  +  k/n) . 

\)-X»  V-K»  0  G 


The  conclusions  of  the  theorem  hold  under  the  more  general  condition  that 
the  distribution  of  the  observations  be  continuous  with  finite  second  moments 


and  support  (-«,<»). 


PROOF  OF  THEOREM  2.  Denote  the  right-hand  side  of  (Al)  by  at  for  t  = 
0,  1,  2,  ...  ,  and  let 


gU)  -  l  a  ?C  for  Id  <  1.  (A2) 

t-0 


Then  a  2  0  for  all  t  and  'E(N)  *  I  a  . 

t-0  C  t 

For  part  (a),  if  c  <_  0  then  $(ci/u-<5)  <f(-6)  for  all  u  ^  0,  and  a£  <_  $  (-6) 

for  all  t.  Since  0  <  $(-<$)  <  1  for  all  S,  £at  is  absolutely  convergent 

and  E(N)  is  given  by  (3.2).  If  0  <  c  <  i/T  then  by  Tauber's  second  theorem. 


lim  g(?)  »  M  <  a  =M. 

G-*-l”  t-0 

Now,  for  sufficiently  small  ]d>  we  find 
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g(C) 


1 

rOsv) 


®  %(v-2)  -u 

^  tl  G  j 

/  - -  du 

0  1-C<S(- 6+c  ni) 


(A3) 


and  we  shall  show  that  g(c)  is  convergent  for  all  |?|<1.  Using  (see 
Bimbaum,  1942) 


l-$(y)  >  - -  for  y  >  o,  (A4) 

y  +  (y  +4)^ 

where  4> (y)  is  the  standard  normal  density,  we  get 

1  -  $(c>/u-<5)  >  - -  if  M  >  s/c. 

C/11-6+  {(c>/u-(5)2+4}^ 

Hence,  if  6 < 0  and  0  s  ;  $  1, 

2  oo  2 

[e**6  ^7/2?r0sv)]  /  [c,/ir-6+{(c^-6)2+4}1'S]uJ5(v'"2)e'u(1_l5C  ^^dus#  <«,. 

0  1 

Thus,  g(^)  exists  for  all  |^|  s  1,  and  (3.2)  follows  by  letting  5  ■+  1  .  If 

<5  >  0,  the  proof  is  similar.  One  now  splits  the  range  of  integration  in  (A3) 
2  2  2  2 

into  [0,6  /c  ]  and  [6  /c  ,<*>],  and  applies  the  technique  above  to  the  second 
integral . 

For  part  (b) ,  by  Abel's  theorem, 

00  00 
J  a  <  *  =£>lim  g($)  =  Z  a  . 
t=0  ^1"  t-0  c 

Using  (see  Poliak,  1957) 
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g(0  >  [e^2/2ir/2r (%v)  ]  /  [cv£-<5  +  {(c^-<5)2  +  8/ir}Ss]e"u(1'Jac  )-<ScvSlu 

0 

The  right-hand  side  diverges  for  5  =  1  and,  since  g(?)  is  continuous  at 

5=1,  lim  g(4)  =  «.  Consequently,  E(N)  =  £a  =  «  for  all  <5  <  0.  The  proof 
C+l" 

is  similar  when  6  >  0. 

For  part  (c) ,  if  <5<0,  we  find  g(l)  =  00  by  the  method  of  part  (b) ,  and  E(N) 
If  (5  >  0,  we  find 


l  u^(v-2)e-u  “ 

gU)  *  TC^[  q  !-;$(- 5fcvG)  du  +  ^2  l-5$(-<^c/I)  du]s  (1-hO  r 


%(\>-2)  -u 
u  e 


f  u%(V~2)e~Udu  +  ^  [v^-^{(v^T-6)2+4}J5]uJ5(v'2)e"'S’^rdu 

o  2cr(l5V)  0 


s  h2  <  <=° 


and  (3.2)  follows. 


To  prove  the  last  part,  we  note  S  *  c  in  probability,  and  h(S)  • 

[1  -  $(-<5  + S)]^is  a  continuous  function.  Consequently,  h(S)  -*■  h(a)  in 

2 

probability.  It  suffices  now  to  show  that  sup^  E(h  (S))  <  ~  (see  Loeve,  1963, 
Corollary  2  on  p.  164).  Using  the  facts  that  1  -  <!>(x)  S  1  -  <S>(x0)  for  x  <  xQ 
and  1  -  <f(x)  >  <|>(x)(x_'L  -  x  )  for  x  >  xQ,  where  xQ  >  0  is  fixed,  it  can  be 
shown  after  some  algebra  that 

[1  .  $(_<5  +  s)]"2  <;  [l-$(2+j  5 1 )  ]-2  if  S<(2+|<5|+<5)<r/k»£ 

2  2 

<  2it  e2<^  +  S  if  S>(2+|  <5  |+<5)a/kv'n 


\  .  —  H  M  +  — ?- 

2  u  2+ 


+6 


r). 


where 
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Since 

E(eXS  )  -  (1-2 Xa2/v)'v/2  <  (2Aa2+l) X(j2+ls 

2  2 
for  v  >  2Aa  +1,  it  follows  that  sup  E(h  (S))  £  M,  where 

2  2 

M  -  [l-<t(2+|  5  j )  ]-2  +  2re25  (l+2Xo2)Xa  +*s . 
Consequently,  E(N)  -*•  [1-$  (-S+kr'n)  ]_L  as  v  ■*  «. 


PROOF  OF  THEOREM  3.  From  (3-1)  we  have 

P(N  >  t)  -  E($C(-<5  +  —  £)). 

a 

Since  <tCis  a  convex  function  of  $,  Jensen's  inequality  gives 

P(N  >  t)  *  E(4(-«  +  M  /S'))t 

«  P((Xi  -  uQ)/S  *  k)C 

The  theorem  follows  easily  from  this  result.  Note  that  the  argument  would 
apply  to  any  distribution  function  so  that  the  conclusion  of  the  theorem  is 
also  valid  for  non-normal  distributions  with  finite  variance. 


PROOF  OF  THEOREM  4.  We  use  the  basic  expression  (3.2)  where  5=0.  Using 
(A 5)  for  1  -  $(c^u)  we  get 


«H>  i  2$^  7  Cc*C  +  (c2u  +  S.)1*]  u‘l(v-2>  e-“d-^2ldu  . 

2  2 

Since  c  u+8/ir  >  naxic  u,3/ir}  for  all  u  ^  0,  the  right  hand  side  above  yields. 


upon  integration,  the  first  inequality  in  (3.3).  Similarly,  using  (A4)  and 


{c2u+4}li<min(2(cv^u  +  (1/c/u)),  2(ci4rH)],  u  s  0, 


in  (3.2)  we  get  the  second  inequality  in  (3.3). 
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PROOF  OF  THEOREM  5.  We  shall  use  the  basic  expression  (3.2)  where  5  >  0. 
Write 


E(N)  -  Ix  +  I2  , 


*s(v-2)  -u 
u  e  du 


where 

62/c2  ^(v-2)  -u,  “> 

I1=(r(%v)}  /  $(5-c/S)  ’  I2  =  {r(iiv)}  ^2/c2  l-Kc^-or* 

ry  r\ 

Since  $(0)  <  $(6-c/u)  <  $(6)  for  all  0  <  u  s  5  /c  ,  we  find 


{$(5>r1P(X2(v)<262/c2)<  I1  £  2P(x2(v)  <  262/c2)  .  (A8) 

fS  1 

Consider  now  1^.  Using  (A4)  and  the  fact  { (c/u- 6)  +4 Y"  <  c/u  -  6+2  for 
2  2 

u  >  6  /c  one  gets 

2 

.  /2tt  e 2<S  ,  *s(v-2)  .-u(l-Asc2)-6c/u  ,  , 

I2  <  -p(W) —  J 2  2  tc/u  +  e  du*  (A9) 

'  "  **  <5  /c 

Using  (A5)  and  the  fact  { (c/u-5)2+8/ir}^  >  max[c/u  -5,  /8/tt]  for  u  >  52/c2 
one  gets 


.  2  „ 

/Hr  ea  ,  w,  x  4(v-2)  -u(l-Jsc  )-  c/u 

X2  1  TftS-  f62/c 2  «»>u 


du,  (A10) 


where  M(u)  =  2(c/u  -6)  or  c/u  -6+/8/tt. 
Inequalities  (A8)  through  (A10)  lead  to  (3.4). 


PROOF  OF  THEOREM  6.  The  result  follows  by  writing  c  =  /2  in  (A8)  through 
(A10)  and  using  the  duplication  formula  for  gamma  functions. 

PROOF  OF  THEOREM  7. 


P5(N  >  t) 


P.(X.  -  Un  £  kS,  i*l»  2,  ...  ,t) 
o  1  U 

00  t  , 

f  ft  $(-6.  +  —  /n)  fc(s)ds. 

0  i-1  1  c  S 


Ws  T“ 

Since  <t> (~<5^  +  —  *n)  is  decreasing  in  6^  the  result  follows.  Note  that  the 
theorem  holds  under  the  more  general  condition  that  the  observations  are 
from  a  location  parameter  family  with  location  parameter  6. 

PROOF  OF  THEOREM  8.  If  either  E(N)  <  «  or  E(N")  <  «  then  E(N’)  <  ®  since 
N'  =  min{N,N"}.  Part  (a)  holds  since  E(N)  <  =  for  c  <  Jl.  Now 

CD  00 

E(Nf)  =  T7U~\  2  /  [$(-<$  +  cv'u)  -  <K-5-cv/u)  ]Cu  U  du  (All) 

t=0  0 

which  is  symmetric  about  6*0.  Moreover,  E(N')  decreased  in  |g|  so  that 
E(N')  is  maximized  when  6=0.  Using  methods  similar  to  the  proof  for  Theorem 
2(b),  part  (b)  can  be  proved  by  showing  E(N')  =  ®  when  c  >  ^2  and  6=0.  When 
6^0  part  (c)  follows  from  the  fact  one  of  E(N)  and  E(N")  i3  finite.  If 
6  =  0  we  must  again  use  methods  similar  to  the  proof  for  Theorem  2.  If  E(N') 

<  °°  then  part  (d)  holds  since  the  summation  and  integration  can  be  interchanged 
in  (All) .  The  proof  for  part  (e)  uses  the  same  methods  as  the  proof  of  Theorem 
7(e). 

PROOF  OF  THEOREM  9.  Following  the  proof  of  Theorem  7, 

00  t  _ 

P  (N'  >  t)  -  /  n  P  ( jx  -  y  I  5  ks | S=s)  f  (s)ds 
0  i*l  i 

Since  P^  ( [ Xi  -  £  ks|s»s)  is  decreasing  in  ,  the  result  follows.  Note 

that  the  theorem  holds  for  any  underlying  distribution  which  is  symmetric 
and  unimodal  since  this  condition  is  sufficient  to  guarantee  that  the  distri¬ 
bution  of  is  symmetric  and  unimodal  and  this,  in  turn,  implies  that 
Pg  (|Xi  -  PqI  s  ks]s»s)  is  decreasing  in  16^. 
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T 


a 

v=29 

v=»69 

V=B00 

0 

2100.0 

1088.6 

740.8 

0.5 

75.9 

53.8 

44.0 

1 

8.0 

6.9 

6.3 

2 

1.2 

1.2 

1.2 

3 

1.0 

1.0 

1.0 

Table  1.  ARL  values  when  a  is  estimated  with  v  degrees  of  freedom, 


U-Uq 

a 

t 

v=9 

v*29 

v=*69 

\)SBOO 

0 

100 

.6900 

.7962 

.8384 

.8736 

0.5 

100 

.2329 

.1774 

.1427 

.1001 

1 

20 

.1375 

.0796 

.0543 

.0316 

2 

2 

.0639 

.0382 

.0307 

.0252 

3 

1 

.0064 

.0025 

.0018 

.0013 

Table  2.  Value  of  P(N>t)  when  a  is 


estimated  with  v  degrees  of  freedom 


SECURITY  CLASSIFICATION  Of  THIS  PACE  fWlibn  Polo  Entered) 

I  REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 


]2.  30 VT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (mid  Subtitle)  _ _ _ 

/*  -  —  — ^ 

.  SHEWHART  X-CHARTS^ITH  ESTIMATED JROCESS 


VARIANCE, 


7.  AUTHOROI  ,  ^  •  /  W 

•  -  ~  i  Jfr/ 

I  C  B.K.  Ghosh,  R./ Reynolds Rui  /  ( 

9.  PERFORMING  ORGANIZATION  NAME  ANO  A00RE5S 

Virginia  Polytechnic  Institute 
Blacksburg,  Virginia  24061 

M.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

U.  S.  Army  Research  Office 
Post  Office  Box  12211 

_ Research  Triangle  Park,  N'C  27709 _ 

14.  MONITORING  AGENCY  NAME  4  AOORESSfli  iron i  Controlling  officii) 


f 


5.  T^PZ  ar  REPORT  6  PERIOO  COVERED 

(  Technical  Repert,  ) 

4.  PERFORMING  PRO,  REPORT  NpB^R 

o71  gz_SD _ 

8.  CONTRACT  OR  GRANT  NUMBERfb.1 


i^yDAAGC29-78-G-^172 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  8  WORK  UNIT  NUMBERS 


IS.  SECURITY  CLASS,  (at  title  report) 


_ Unclassified _ 

isT  declassification/ downgrading 

SCHEOULE 

NA 


114.  DISTRIBUTION  STATEMENT  (at  t> )T«  Repot!) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (at  the  abstract  entered  In  Block  20,  It  dllterenl  tram  Report) _ 

;'i? '  .  /.f7v  /'/ 

NA  -  .  -  .  V  7  -  '  ' 

•  V— —  — .  ■•••>• 

t®.  Supplementary  notes 

The  findings  in  this  report  are  not  to  be  c-'r,s  trued  as  an  official 
Department  of  the  Army  position,  unless  so  iesignatei  by  other  authorized 
_ documents .  _ _ _ _ 

19.  KEY  WOROS  (Continue  on  reeeree  aide  it  neceeemry  and  Identity  by  block  number) 

X-Chart;  Shewhart  Chart;  Process  Control;  Control  Chart;  Estimated 
Variance;  Economic  Model;  Average  Run  Length. 


P 5T  ABSTRACT  (ConUmn  on  rororto  cldo  It  neeeeeery  and  Idanilty  by  block  number) 


Properties  of  the  Shewhart  X-chart, for  controlling  the  mean  of  a 
process  with  a  normal  distribution  are  investigated  for  the  situation 
where  the  process  variance  must  be  estimated  from  initial  sample 
data.  The  control  limits  of  the  X-chart, depend  on  the  estimate  of  d4 

2  /  .  r  *  I 

and  thus,  unlike  the  case  when  a  is  known,  the  X-chartA is  not  equivalent 


DO  1473  EOITION  OF  *  NOV  «8  IS  OBSOLETE 


Unclassified  ^  A-4-/" 

SECUBITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  inland) 


/  Jr 


y  o  7  s  u 


SECURITY  CLASSIFICATION  OF  THIS  PACEfTWin  Data  Bnfnd) 


ZU.  , 

-  £  -  'jua<'T»< 

2" 

to  a  sequence  of  independent  tests-  When  a  is  estimated  the  distribution 

of  the  run  length  is  not  geometric  and  cannot  be  characterized  simply  in 

terms  of  the  probability  of  a  signal  at  a  given  point.  The  average  run 

length  (ARL)  for  the  X-chart,,  is  expressed  in  terms  of  an  integral  involving 

the  normal  cdf,  and  it  is  shown  that  the  chart  signals  with  probability  one, 

2 

but  the  ARL  may  not  be  finite  if  the  size  of  the  sample  used  to  estimate  c 

is  sufficiently  small.  In  addition,  certain  bounds  for  the  ARL  are  also 

derived.  Numerical  integration  is  used  to  show  that  the  effect  of  using 

2 

small  sample  sizes  in  estimating  o  is  to  increase  the  ARL  and  the  vari¬ 
ance  of  Che  run-length  distribution. 


Jl  URITY  Cl  »SSIFICAT  On  OF  This  PAGEflFhen  Oat*  Enter*)*; 


